Weak ferromagnetism in antiferromagnets: Fe 2 C>3 and La 2 Cu04 
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The problem of weak ferromagnetism in antiferromagnets due to canting of magnetic moments 
was treated using Green's function technique. At first the eigenvalues and eigenfunctions of the 
electronic Hamiltonian corresponding to collinear magnetic configuration are calculated which are 
then used to determine first and second variations of the total energy as a function of the magnetic 
moments canting angle. Spin-orbit coupling is taken into account via perturbation theory. The 
results of calculations are used to determine an effective spin Hamiltonian. This Hamiltonian can 
be mapped on conventional spin Hamiltonian that allows to determine parameters of isotropic and 
anisotropic (Dzyaloshinskii-Moriya) exchange interactions. The method was applied to the typical 
antiferromagnets with weak ferromagnetism Fe2C>3 and La2Cu04. The obtained directions and 
values of the magnetic moments canting angles are in a reasonable agreement with experimental 
data. 
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Heisenberg Hamiltonian is a basis of most theo- 
retical investigations of transition metal compounds 
magnetismii 2 ^ The essential part of these investigation 
is determination of exchange interaction parameters Jy . 
It can be done in a phenomenological way by fitting those 
parameters to reproduce experimental data (temperature 
dependence of magnetic susceptibility and magnon dis- 
persion curves obtained in inelastic neutron scattering 
measurements) i* However much more physically appeal- 
ing is to obtain them in ab-initio calculations. In most 
cases it was done via calculating total energy values for 
various different magnetic moments configurations. Map- 
ping on Heisenberg Hamiltonian gave a system of linear 
equation for Jjj (for example see Ref. |5() . This procedure 
becomes inconvenient for the systems with a large num- 
ber of long range competing exchange interactions like in 
(VO) 2 P 2 7 , NaV 2 5 , Cu 2 Te 2 5 X 2 (X=Br,Cl), etc& 

In 1987 A.I. Lichtenstein et ali proposed the cal- 
culation method that does not use total energy differ- 
ences. They determined exchange interaction parameters 
via calculating second variation of total energy S 2 E for 
small deviation of magnetic moments from the collinear 
magnetic configuration. The expression for this second 
variation was derived analytically and required for its 
evaluation calculation of the integral over production of 
one-electron Green functions. This method was since 
then successfully applied to the various transition metal 
compounds 

The combination of low symmetry and spin-orbit cou- 
pling was shown by Dzyaloshinski: 12 and Moriya™ to 
give rise to anisotropic exchange coupling. Moriya has 
shown how the processes involving an additional vir- 
tual transition due to spin-orbit coupling can cause an 
anisotropic exchange interaction as a correction to the 
isotropic Anderson superexchange term and introduced 
new term in spin Hamiltonian which is Dzyaloshinskii- 
Moriya interaction (DM). Solovyev et ali^ had shown 
that Dzyaloshinskii-Moriya interaction parameters can 



be calculated using perturbation theory and Green's 
function technique and described the canting of mag- 
netic moments of LaMn.03. Recently, Katsnelson 
and Lichtensteini^ derived the general expression for 
Dzyaloshinskii-Moriya interaction term in LDA++ ap- 
proach. 

This paper is devoted to the problem of first-principles 
theoretical description of weak ferromagnetism in antifer- 
romagnets, specifically to the task of calculating weak fer- 
romagnetic moment value and direction of spin canting. 
For this we consider first and second variations of total 
energy of the system at small deviation of magnetic mo- 
ments from collinear configuration with spin-orbit cou- 
pling introduced as a perturbation using Green's function 
technique. We show that there is an additional on-site 
term that was not taken into account in previous work^4 
which gives significant contribution to weak ferromag- 
netic moment. Basing on the results of the calculations 
we propose effective single site Hamiltonian. This Hamil- 
tonian is sufficient for solving the problem of spin cant- 
ing but it also can be rewritten to the conventional form 
containing isotropic and anisotropic exchange interaction 
terms. We have applied our method for weak ferromag- 
netism in Fe 2 03, the classical system which was used by 
Moriya in his pioneering work^ and in antiferromag- 
netic cuprate La 2 Cu04 in low temperature orthorhom- 
bic phase, estimated ferromagnetic moments values on 
the metallic ions in these compounds and determined the 
plane of spin canting. 

Briefly, this paper is organized as follows. In Sec. I 
we describe the method for calculation spin Hamiltonian 
parameters responsible for magnetic moments canting. 
Sec. II contains the results of our calculations for Fe 2 03 
and La 2 Cu04 crystals. In Sec. Ill we discuss and briefly 
summarize our results. 
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I. METHOD 

According to the Andersen's "local force 
theorem"i§iiL±& the total energy variation SE un- 
der the small perturbation from the ground state 
coincides with the sum of one-particle energy changes for 
the occupied states at the fixed ground state potential. 
In the first order for the perturbations of the charge and 
spin densities one can find the following relation^: 



/Ef rEp 
de e 5n(e) = ep Sz — / de 
-co J — oo 



5N{e) 



de6N(e), (1) 



here n(e) = dN/de is the density of electron state, 
N(e) is integrated density of electron state and ep is 
Fermi energy. In the case of magnetic excitation the 
change of total number of electrons Sz equals zero. The 
Green function G is formally expressed in the usual way 
G = (e — One can express density of states and 

integrated density of states via Green function G: 



n(e) 



Ira Sp G(e) 



(2) 



and 



New Hamiltonian of the system after rotation of the spin 
on j site around direction n on the angle \8cj>\ 



H = WHU 



(9) 



The first variation over the angle of rotation is expressed 
in following form: 



SH = - i 60 [H, d] 



(10) 



In the basis \ilma) (where i is site, I is orbital quan- 
tum number, m is magnetic quantum number and a 
is spin index) the Hamiltonian matrix takes the form 
Hum jim' = (ilTrw\H\jlm'a J }. For simplicity below we 
drop the index of orbital and magnetic quantum numbers 
and leave spin and site indexes. We assume that with- 
out spin-orbit interaction the ground state corresponds 
to the collinear magnetic configuration at which all spin 
moments lie along z axis. Therefore the Hamiltonian ma- 
is diagonal in the spin subspace 



trix Hf/ 




Hi 



One can rewrite the first variation of Hamiltonian Eq. i|10|) 
in the following form 



N(e) 



Ira Sp ln(e — H). 



(3) 



Then the variation of integrated density of state is given 
by 



1 



SN(e) = -ImSp[SHG]. 



(4) 



Therefore the first and second variations of total energy 
of the system take the following forms: 



SE = - 1 - 

7T 



Ef 



deImSp(5HG) 



and 



5 E = 



de Ira Sp (5 HG + SH G SH G). (6) 



Operator of spin rotation on the site j on the angle 



1 5(f> | around direction n = —$- is given by 



U = 



(7) 



where a = (a x , a y , a z ) are Pauli matrices. For small 

| S(j) | values we can expand the spin rotation operator in 
following way 



U 



1 



1 



1 



(8) 



SHi 



( 
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A, \ 


1 Aj 


l) 




5) 


V 2 


V 2 



(11) 



where A., = Hi. — H.. . It is easy to show that the 
second variation of Hamiltonian is given by 



S 2 H 3J = S\ 



2 

u 2 





A, 



s2 *u „- ¥ 



2 







(12) 



The rotation of spin moment around z axis doesn't 
change the energy of the system therefore the term with 
(5) 5(j) z is absent in Ea. [|llll2[) . 

Then we take into account the spin-orbit coupling via 
perturbation theory. The Green function in the first or- 
der perturbation theory with respect to the spin-orbit 
coupling can be written as 



Gij — Gij 



Hl° G 



kj - 



(13) 



here H%° = Xf. L S, i, j and k denote site, G 4 -? is Green 
function of system with collinear magnetic configuration, 
and Afc is spin-orbit coupling constant. The first variation 
of total energy Eq.© takes the form: 



SE = 



7T ^ ' 



E F 



de Ira Sp 



: {SHi Gu + ^ SHi Gik H^ Gki)- 



(14) 
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The first term in Eg. i|14|) is zero. The second term can 
be expressed as a following sum 



SE = Y^ A i^i +^6$, 



(15) 



where 



E F 



Y. ( a < <4 ffJI , Gl - A, G]„ H- t ° n G s kl ) (16) 



and 



/to 



2 (Ai Gj fe tff iT + Aj <4 fff n Gij. (17) 



We consider the situation when all spins lie along z axis 
and therefore the rotation around it does not change 
the energy of system. In order to find A| component 
of the magnetic torque vector Ai we change the coordi- 
nate system in the following way (x,y,z)— >(z,y,-x) (rota- 
tion around y axis): 

--K-.i)(2S5§)(-)- <■« 

Therefore Af component in new coordinate system is Af 
in the old one: 



i k 



4tt 



de Re 



X ( Ai G *k ( H k°n ~ H k°lO G L 

k 



(19) 



In contrast to the first variation SE, the second varia- 
tion of total energy S 2 E for small deviations of magnetic 
moments from ground-state collinear magnetic configura- 
tion has nonzero value without taking into account spin- 
orbit coupling: 

1 P Ef 1 
5 2 E = / deImSp(- ^ViT^ 

+ lH PHiiGii + E SH * (: :> 6,1 ■ (: - ( 2 °) 



where 



Sp(8 2 H vl G u ) = -5 2 0f Ai(G^ - 



Is^MGi - Gl 



(21) 



and 



S P (6Hi G {j SHj Gft) = - 8tf 8$ (A, Gj, A, G T 



-8<%5<j>) 



ij ^3 ji 

1 a . n\ 



(A i Gl i A j Gj i ).(22) 



Using the condition G\ - G l u = (G T AG^ = 
J2j ^JjAj Gjj one may rewrite Ea. H2U|) in following form: 



S 2 E = — 
4ir 



d< Im ^(A,G;,A,G;, 



x({Stf -S^Y + (8$ -84>]Y). (23) 

One can see that the Eq. (|23|l contains only x and y com- 
ponents of 8<j). In order to include z component one can 
use the same rotation of coordinate system as for the site 
magnetic torque vector A Eq. (|18|) . Finally, we obtain the 
following function of the total energy over angle 8<p: 



AE = J2 A * S & + lj2 J y 



M 2 , 



where 



Ji i — 



(24) 



(25) 



The aim of this paper is description of canted mag- 
netism in transition metal compounds caused by spin- 
orbit coupling. For this we have used the expression 
Ea. (|24|) for the total energy as a function of canting an- 
gle. In order to solve the problem of the weak ferromag- 
netism in antiferromagnets we suppose that the crystal 
is an antiferromagnet containing two sublattices 1 and 2, 
with the same canting angle for the atoms belonging to 
the same sublattice. With this assumption the Eq.J53J 
is reduced to the following form: 



AE = A 1 



A 9 



8fo\ 2 . (26) 



Our results for Fe2 03 and La2Cu04 demonstrated that 
A\ = - Ai (torque vector has an opposite sign for the 
atoms belonging to the different sublattices). That gives: 



AE = A x 



(27) 



If we further suppose that the deviations of magnetic mo- 
ments from the average direction defined by 8<pi have the 
same absolute value but different sign for both sublat- 
tices, then Ea. (|27l) takes the following form (we suppose 
that magnetic moments lie in plane perpendicular to site 
magnetic torque vector A and canting occurs in the same 
plane) 



AE = 2Axi 



h + 4 X J u Nil 2 - 



(28) 
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Then we find the value of \S(p±\ where AE has a mini- 
mum: 



\84>i 



(29) 



The next step is to establish a connection between 
Eg.l|24p and conventional spin Hamiltonian 

H = Hdm + Hexch = E Dij [Si x ej] 

+ E ■'•'''■''>■ ( 3 °) 

i¥=3 



where ei is unit vector in the direction of the ith site 
magnetization, is exchange interaction and is 
Dzyaloshinskii-Moriya vector. One can rewrite the sec- 



ond term in Ea. (|3U[l as H e 



e, cos( 



In the limit of small canting angle values we can assume 



that 



= 1- \Hi 



2 /2 and exchange inter- 



action energy for antiferromagnetic configuration has a 
form: 



AH, 



(31) 



Therefore in the limit of small 5<j) we can directly map 
the second term of total energy variation Ea. (|24|l onto 
first term in spin Hamiltonian Eci. (|3U[l . 

The first term in Ea. H24|) describes the deviation of 
the spin moment on the site i from the initial collinear 
spin configuration direction. We assume that this initial 
spin direction on the site i is defined by the direction of 
Weiss mean-field Hf rF — J2j(^i) J*j ^3 (corresponding 

-). Therefore we can map the 



unit vector is e„ 



Hi 



first term in Eq.JHJ on the spin Hamiltonian 



H, 



dec 



£4 



(32) 



describing the deviation of spin moments away from the 
direction of external Weiss field. (We have used 



here the connection between rotation vector 



and the 



change of the magnetic moment unit vector Set = €i — e®\ 
S&i = [64>i x e^ ].) In order to demonstrate the connection 
between A, and one can rewrite the Ea. (|32|l in the 
following form: 



H de 



1 



2 VZ^ 

i 

Y, 1 



A; 



ttWF 

i 

WF\ 



\H 



H WF 
\Hf F \ 



x ej]). 



(33) 



Using our definition of H^ F we obtain: 

Hdev = \(Y.T^F-M^ X 



\3WP\ 
ij I I 



\H WF 

y 1 3 



(34) 



This gives us the following expression for parameter 
of spin Hamiltonian Ea. l|3U|) : 



2 J%3 ^\H WF \ 



A, 



\H 



WF 



0- 



(35) 



Therefore the components of Dzyaloshinskii-Moriya in- 
teraction vector are given by 



D x = -— J 



deRe ^ 



X( 



kj 



\Hf F \ 



Ai G.j. Hf° n GL — AjGJj. -fff ti G 



7 ik k TT hi 



-4 



| J. fci • 



D y = — 7 



(36) 



de Im 



k 



^3 G jk H k°n G lj + AjG^ H k°n G l 



kj 



\Hf F \ 

^■j G ik H k°n G L + ^i G lk H k°n G L 
\hF f \ 



(37) 



A, G ]k 


( TTSO 




Gl k3 










Aj G jk 










\Hf F 






^ G l 


( TTSO 

sP-k TT 


rrso ' 

U kU- 


) G l 



+ 



\ H WF\ 
Ai G \ k (Hk TT ~~ ^lii) 



; TT 

\H 



WF\ 



(38) 



We have obtained more general expression for 
Dzyaloshinskii-Moriya interaction parameter in compar- 
ison with those presented in paper There are two kind 
of contributions into magnetic torque vector Af. on-site 
interaction Bn (absent in workA^) and intersite interac- 
tion Bik (i ^ k) . We have found that on-site contribution 
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in magnetic torque A which was not considered before 
plays an important role in weak ferromagnetism descrip- 
tion. 

We have applied the calculation scheme developed 
above to the typical antiferromagnets with weak ferro- 
magnetism Fe2C>3 and La2CuC>4 in low-temperature or- 
thorhombic phase. In order to calculate Green functions 
corresponding to the collinear spin configurations we used 
LDA+U approaches realized in LMTO method within 
Atomic Sphere Approximation^ 



II. RESULTS 
A. Fe 2 3 

Weak ferromagnetism or weak non-collinearity of es- 
sentially antiparallel magnetic moments was first ob- 
served in ev-hematite, a-Fe203. 2 « The trigonal crystal of 
Fe2C>3 has R3c space group. Depending on temperature 
a-hematite may occur in two different antiferromagnetic 
states: at T < 250 K the spins are along the trigonal axis, 
and at 250 K < T < 950 K they lie in one of the verti- 
cal planes of symmetry making a small angle 1.1 x 10~ 3 
with basal planei 2 ^* 2 ^ In the latter case the a-Fe2 03 has 
a net ferromagnetic moment. Dzyaloshinskii has shown 
that the spin superstructure gives rise to a nonvanish- 
ing antisymmetric spin coupling vector which is paral- 
lel to the trigonal axis. Moriyai^ gave phenomenologi- 
cal Dzyaloshinski's explanation a microscopic footing by 
means of Anderson's perturbation approach to magnetic 
superexchange. 

Sandratskii et alS^ have performed the calculation 
based on the local approximation to spin-density func- 
tional theory (LSDA) using the fully relativistic version 
of ASW method. In spite of the well-known problem 
that the LSDA has with a proper determination of the 
energy gap in semiconducting and insulating materials, 
the authors 2 * succeeded in describing a weak ferromag- 
netism and obtained the ferromagnetic moment of about 
0.002 /i^. In the present study we treat the problem of 
description of weak ferromagnetism in Fe2C>3 using per- 
turbation theory. 

Electronic structure of a-hematite calculated using the 
standard LDA+U approximation 1 - with on-site Coulomb 
interaction parameters U = 5 eV, J = 0.88 eV and struc- 
ture data fromS^ is in a good agreement with previous 
theoretical calculations^ We have obtained the mag- 
netic moment of 4.1 hb per Fe atom. This value is a 
little smaller than those obtained in experiments (4.6-4.9 
/i_b). The energy gap value of 1.67 eV is also slightly 
underestimated comparing with experimental data (2.14 
eV in paper—). 

The Brillouin zone integration has been performed in 
the grid generated by using (6;6;6) divisions. The energy 
integration has been performed in the complex plane by 
using 800 energy points on the rectangular contour. The 
calculated isotropic exchange interactions and contribu- 



^ * X 
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FIG. 1: Crystal structure of a-Fe203. Large circles are Fe 
atoms which belong the smallest unit cell used in the LDA+U 
calculations. The small circles are Fe atoms which surround 
Atom 2. Arrows denote the magnetic configuration used in 
our calculation. 

tions in site magnetic torque A 2 are presented in Table I. 
The simplified crystal structure and the interaction paths 
are shown in Fig.l. 

One can see that the obtained interaction picture is 
more complicated that those Moriya examined in order to 
describe the weak ferromagnetism in a-hematitei 3 There 
are strong isotropic exchange interaction with third and 
fourth neighbours. This agrees well with experimental 
results 28 and theoretical predictions^ The summary ex- 
change of atom 2 with nearest neighbours is given by 
^2 = £^ 2 J 2J = 189.26 meV. 

The components of site magnetic torque vector on the 
atom 2 are given by A% — £\ B 2i = e V, A\ — 
£\ B v 2i = eV and A\ = Y.i B li = 0.281 meV. One can 
see that the on-site interaction _B| 2 gives the main contri- 
bution in A2. We have calculated also the site magnetic 
torque of A\ which has the following components (0;0;- 
0.281), the same value but the opposite sign comparing 
with A 2 . The value of canting angle of \6<p\ — 0.4 x 10 -3 
calculated with Ea. (|29|l is of the correct order of magni- 
tude but is more than two time smaller than experimental 
data 1.1 x 10~ 3 fRef. l2l2^ . The reason for the difference 
between experimentally observed and calculated here val- 
ues of canting angle could be the necessity to take into 
consideration the higher order terms with respect to spin- 
orbit coupling which were not considered in the present 
study. 

It is easy to show that in the case when all spin lie 
along z axis there is no canting of the spin moments. On 
the other hand if direction of Weiss field is perpendicular 
to z axis the canting exists and the system has weak 
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TABLE I: Isotropic exchange interaction parameter and dif- 
ferent contributions in components of site magnetic torque of 
Fe2C>3 (in meV). The couplings with negative sign are ferro- 
magnetic, dij is distance between i and j atoms in a.u. Rij 
is radius vector from i site to j site in units of the lattice 
constant (5.49 a.u.). 





dij 


Rij 










(2,2) 





(0;0;0) 











0.162 


(2,1) 


5.45 


(0;0;-0.99) 


8.576 








0.005 


(2,3') 


5.60 


(-0.5;-0.86;0.20) 


-7.3 


-0.036 


0.015 


0.001 


(2,3') 


5.60 


(1;0;0.20) 


-7.3 


0.032 


0.023 


0.001 


(2,3') 


5.60 


(-0.5;0.86;0.20) 


-7.3 


0.004 


-0.038 


0.001 


(2,1') 


6.36 


(0.5;-0.86;-0.58) 


25.224 


0.071 


0.019 


-0.14 


(2,1') 


6.36 


(-l;0;-0.58) 


25.224 


-0.052 


0.052 


-0.14 


(2,r) 


6.36 


(0.5; 0.86;-0.58) 


25.224 


-0.019 


-0.071 


-0.14 


(2,4') 


6.99 


(0.5;-0.86;0.79) 


17.502 


0.168 


0.063 


0.101 


(2,4') 


6.99 


(-1;0;0.79) 


17.502 


-0.139 


0.114 


0.101 


(2,4') 


6.99 


(0.5;0.86;0.79) 


17.502 


-0.029 


-0.178 


0.101 


(2,4") 


6.99 (-0.5;-0.86;-0.79) 


17.502 


0.128 


0.094 


0.076 


(2,4") 


6.99 


(l;0;-0.79) 


17.502 


0.017 


-0.158 


0.076 


(2,4") 


6.99 


(-0.5;0.86;-0.79) 


17.502 


-0.145 


0.064 


0.076 



ferromagnetic moment. This picture fully agrees with 
experimental and theoretical datapi 3 * 2 ^ 3 " 



B. La 2 Cu0 4 

In the case of the cuprates Dzyaloshinskii-Moriya in- 
teraction is the leading source of anisotropy, since single- 
ion anisotropy does not occur due to the S=^ nature of 
the spins on the Cu 2+ sites. The experimental data2S*2i 
demonstrate that in case of low temperature orthorhom- 
bic phase the spins do not lie exactly in the Cu-0 
planes, but are canted out of the plane by a small an- 
gle (0.17°). Coffey and coworkers 32 made complete ex- 
amination of the anisotropic exchange interaction in or- 
thorhombic phase based on a symmetry consideration. 
They assumed rotation axis of the the CuOe as a direc- 
tion of antisymmetric exchange interaction and obtained 
that the spins are canted in plane which is perpendicular 
Dzyaloshinskii-Moriya vector. 

The first attempt at a microscopic calculation of the 
Dzyaloshinskii-Moriya anisotropy for La 2 Cu04 in low 
temperature orthorhombic and tetragonal phases was 
made by Coffey, Rice, and Zhangi^ 3 . They have neglected 
the symmetric anisotropic exchange interaction which is 
of the second order with respect to the spin-orbit coupling 
and can be written in spin Hamiltonian in the following 
form: 



Ha — Si 



where My- is a symmetric 3x3 tensor. 

Shckhtman, Entin-Wolhman, and AharonySi have 
shown that symmetric anisotropic exchange interaction 
can not be neglected, since its contribution to the mag- 
netic energy is of the same order of magnitude as that 



f t 

O O 

4 3 

t t 



FIG. 2: The schematic crystal structure of La2Cu04 in low 
temperature orthorhombic phase. The open circles are oxygen 
atoms which are tilted up out of the Cu plane, the black circles 
are oxygen atoms tilted down out the Cu plane. The big 
gray circles are copper atoms. Arrows denote the magnetic 
configuration used in LDA+U calculations with spin moments 
lie along z axis. 



of the antisymmetric anisotropic exchange interaction 
(Dzyaloshinskii-Moriya). On the basis of our results (as 
we show below) we can conclude also that taking into 
account the second order terms with respect to the spin- 
orbit coupling must be important. 

We have performed the LDA+U calculations for 
La2Cu04 in low temperature orthorhombic phase using 
structural data for Nd doped La2CuC , 4j 3 £ with on-site 
Coulomb interaction parameters U = 10 eV, J =1 eV 
(the same as used in workSi). The schematic structure of 
Cu-0 layer of La2Cu04 in low temperature orthorhombic 
phase is presented in Fig. 2. 

The experimental value of the energy gap is reported to 
be about 2 eV (Ref.f3?j). Our gap value of 1.94 eV is in a 
good agreement with experimental data. The calculated 
magnetic moment value on Cu atom is 0.61 fiB which 
also agrees well with experiment^ 



TABLE II: Isotropic exchange Joj and the components of dif- 
ferent contributions in site magnetic torque Ao (in meV). Rij 
is radius vector from i site to j site in units of the lattice 
constant (10.14 a.u.). 



(i, ,i) 


Rij J ij 


Dij 


B H 

« 


B tj 


(0,0) 


(0;0;0) 


0.101 








(0,1) 


(-0.5;0.5;0) 14.576 


0.020 


-0.032 


-0.005 


(0,2) 


(0.5;0.5;0) 14.576 


0.020 


0.032 


0.005 


(0,3) 


(0.5;-0.5;0) 14.576 


0.020 


-0.032 


-0.005 


(0,4) (-0.5;-0.5;0) 14.576 


0.020 


0.032 


0.005 



We have performed calculations of isotropic exchange 
interactions and the different contributions to site mag- 
netic torque components (Table II) using the energy in- 
tegration in the complex plane with 700 energy points on 
the rectangular contour and the Brillouin zone integra- 
tion has been performed in the grid generated by using 
(6;6;6) divisions. The obtained values of exchange inter- 
action parameters are in a good agreement with results 
of previous calculations for low temperature tetragonal 
phased and experimental estimations^ The exchange 
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interactions with next neighbours are negligibly small. 
The summary exchange and the components of site mag- 
netic torque are given by Jo = J2i^oJoi — 58.304 meV, 
Al = 0.18 meV, A% = meV and A% = eV. We 
obtained that A2=(-0.18;0;0), again of the same value 
but the opposite sign comparing with Aq. It means 
that the system has net ferromagnetic moment if spins 
lie in plane which is perpendicular to x axis, which is 
axis of rotation of oxygen octahedra. This fully agrees 
with results of previous theoretical works. 30 32-33 The ob- 
tained value of canting angle \5<f>\ = 0.77 x 10~ 3 is about 
three times smaller t han th ose experimentally observed 
2.2 -r 2.9 x 10~ 3 (Ref- HEH) . 

Again as for Fe2C>3 we expect that the inclusion of 
second order spin-orbit coupling terms could improve the 
agreement. Such calculations are in progress. 

III. CONCLUSION 

We present a method for calculation of spin Hamilto- 
nian parameters responsible for magnetic moments cant- 



ing. The effective Hamiltonian for canted magnetism 
was proposed. We shown that the parameters of this 
model Hamiltonian can be obtained from first-principles 
calculations. Using the developed method we described 
the weak ferromagnetism in Fe2C>3 and La2Cu04. It 
was shown that on-site contribution Ba in site magnetic 
torque Ai plays the crucial role for net ferromagnetic 
moment of Fe203 and La2Cu04 in low temperature or- 
thorhombic phase. 
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